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Magnetization hysteresis loops and the ac susceptibiUty, X = x' + ix"i of ^ superconducting thin 
disk are calculated in the critical-state model assuming a field-dependent critical current density, 
Jc{B). The results are obtained by solving numerically the set of coupled integral equations for the 
flux and current distributions [PRB 60, 13112 (1999)] for a disk placed in a perpendicular applied 
field Ba- From the magnetization curves the range of fields where the vertical width of the loop, 
AAI{Ba), relates directly to Jc{Ba) is determined. The susceptibility is analyzed in the limits of 
small and large ac-field amplitudes Bam, and also as a parametric relation x"{x')- Comparing our 
results with experimental data for x"(x') shows that by taking the B-dependence of Jc into account 
the agreement improves dramatically, in particular at small |x'| (large field amplitudes). We show 
that the asymptotic behavior for large Bam changes from x' k Bam and x" o<; Bam for the Bean 
model, to x' oc Bam and x" oc Bam for Jc decreasing with \B\ as \B\~^ or faster. For small Bam the 
behavior can always be described by an effective Bean model with a renormalized Jc- We also find 
that in the x"(x') pfo^ the peak of x" increases in magnitude and shifts towards x' = when Jc 
decreases with \B\. This allows an easy experimental discrimination between a Bean model behavior, 
one with Jc{B), and one where flux creep is an ingredient. 
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I. INTRODUCTION 

The critical state model (GSM) is widely accepted as 
a powerful tool in the analysis of magnetic properties of 
type-II superconductors. In the parallel geometry, i.e., 
for long samples like slabs and cylinders placed in a par- 
allel magnetic field, an extensive amount of theoretical 
work has already been carried out.|-|Exact results for fliis 
density profiles, magnetization,trcl ac susceptibilityETEl 
etc., have been obtained for a number of different field- 
dependent critical current densities. During the last 
years even more attention has been paid to the CSM anal- 
ysis in the perpendicular geometry, i.e., for thin samples 
in perpendicular magnetic fields. Assuming a constant 
critical current (the Bean model), explicit ajjialytical re- 
sults have beenjobtained for a long thin stripLfQ and a thin 
circular disk.Erll3 From experiments, however, it is well 
known that also in such samples the critical current den- 
sity jc usually depends strongly on the local flux density 
B. Due to the lack of a proper theory, this dependence 
often hincLera_a precise interpretation of the measured 
quantitiesll3il3 

In the perpendicular geometry, the ac susceptibility be- 
yond the Bean model has beep-calculated only by carry- 
ing out flux creep simulationgl3'E£l assuming a power-law 
current- voltage relation with a large exponent. However, 
quite recently an exact analytical approach was_devel- 
oped for the CSM, analysis of a long thin stripc3 and 
thin circular disk.E2l In both cases a set of coupled inte- 
gral equations was derived for the flux and current dis- 
tributions. In the present paper we solve these equa- 
tions numerically for the thin disk case, and calculate 



magnetization hysteresis loops as well as the complex ac 
susceptibility. Results for several commonly used jc{B) 
dependences are presented. 

The paper is organized as follows. In Sec. O we give a 
short description of the exact solution for the disk prob- 
lem. In Sec. Ill, magnetization hysteresis loops are cal- 
culated and the relation between the width of the loop 
and jc is discussed. The results for the complex ac sus- 
ceptibility are presented in Sec. IV and analysed with 
emphasis on the asymptotic behavior at small and large 
field amplitudes. Finally, Sec. V gives the conclusions. 



II. EXACT SOLUTION 



Consider a thin superconducting disk of radius R and 
thickness d, where d <^ R. We assume either that d > A, 
where A is the London penetration depth, or, if d < A, 
that A^/d ^ R. In the latter case the quant ityj—X^/d 
plays a role of two-dimensional penetration depthEj We 
put the origin of the coordinates at the disk center and 
direct the z-axis perpendicularly to the disk plane. The 
external magnetic field B^ is applied along the z-axis, 
and the z-component of the field in the plane z = is de- 
noted as B. The current flows in the azimuthal direction, 
with a sheet current denoted as J{r) = J_^,^j{r,z)dz, 
where j is the current density. 



A. Increasing field 

We begin with a situation where the external field Ba 
is applied to a zero-field-cooled disk. The disk then con- 
sists of an inner flux-free region, r < a , and of an outer 
region, a < r < R, penetrated by magnetic flux. 

In the CSM with a general Jc{B) the current and flux 
density distributions in a disk are given by the following 
coupled equationsEil 



J[r) 



-^ / dr' 



-Jc[B{r) 



-r2 MB{r')] 



r < a 



a <r < R 



B{r) = Ba 



Mo 

27r 



F{r,r')J{r')dr' . 



Ba^ 



A*o 



dr 



Vr 



-:Jc[B{r)]. 



(1) 



(2) 



(3) 



Here F{r,r') ^ K{k)/{r + r') - E{k)/{r - r'), where 
k{r,r') = 2vrr^/(r -|- r'), while K and E are complete 
elliptic integrals. In the case of constant Jc, these equa- 
tions reduce to the exact Bean-model formulas derived 
in Refs. [lO| and |ri|. 

Note that the calculation can be significantly simpli- 
fied at large external field where a — > 0, and the critical 
state J{r) = Jc[B{r)] is established throughout the disk. 
The distribution B{r) is then determined by the single 
equation 



B{r) =Ba-^ j Fir, r')J,[B{r')]dr' , 



(4) 



following from Eq. 



B. Subsequent field descent 

If Ba is reduced, after being first raised to some max- 
imum value Bann the flux density will decrease in the 
outer part, a < r < R, and remain trapped in the inner 
part, see Fig. ||. We denote the flux front position, the 
current density and the field distribution at the maxi- 
mum field as Om, Jmir) and Bm{r), respectively. Evi- 
dently, Jm(r), Bjn{r), and a^ satisfy Eqs. (Q)-©. 

Let the field and current distributions during field de- 
scent be written as 

B{r) = B„,{r) + B{r), J{r) = J,,,{r) + J{r). (5) 
The relation between B{r) and J(r) then readsEil 



Jiv) 



TT 



Jc{r), 



la^- 


-r2 J,{r') 




r'2. 


- a2 r'2 - r2 ' 


r < a 
a <r < R 



B{r) = Ba 



Ba 



Mo 

27r 



(6) 



F{r,r')J{r')dr' . (7) 



where we defined 



Mo /-^ Jc{r) 
2 J a ^r^^^ 



dr , 



Jc{r) = Jc[B„^{r) + Bir)] + Jc[B„,{r)] 



(8) 



(9) 



Again, setting J^ ^cari^, these equations reproduce 
the Bean-model resultsJij'M 

If the field is decreased below —Bam the memory of 
the state at Bam is completely erased, and the solution 
becomes equivalent to the virgin penetration case. If the 
difference Bam ~ Ba is large enough then one can again 
use Eq. (Q), only with the opposite sign in front of the 
integral. 

Given the Jc (i?)-dependence, a complete description of 
any magnetic state is now found by solving the equations 
numerically. An efficient iteration procedure is described 
in Ref. ^ 




FIG. 1. Flux density profile as the applied field descends 
from a maximum value B„ „, . 



III. MAGNETIZATION 



The magnetization of a disk is defined as the mag- 
netic moment, tt L r^ J{r) dr, per unit volume. Due 
to symmetry the magnetization is directed along the z- 
axis. In a fully penetrated state described by the Bean 
model with critical current J^q, the magnetization equals 
A/o = JcoR/id. It is convenient to use Mo for normal- 
ization, i.e. 



3 /-^JW 2 



M 



Mo i?3 7o Jc 



r dr. 



JcO 



(10) 



The magnetization can be calculated using the current 
profiles obtained by the procedure described in the previ- 
ous section. Shown in Fig. are magnetization hysteresis 
loops calculated for the Jc (i3)-dependences: 
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FIG. 2. Magnetization hysteresis loops for a thin disk for 
the Bean model {Jc ~ const), the Kim model, Eq. (Ill), 
and the exponential model, Eq. (|l2|), the last two both with 
Bo = 3Bc- The parameters Be and Mo are defined in the 
text. 



Jc - Jco/il + \B\/Bo) 
Jc = JcO exp(-|B|/Bo) 



(Kim model), (11) 

(exponential model). (12) 



A striking manifestation of the i?-dependence is a peak 
occuring at small Ba ■ The calculations show that for any 
choice of the parameter Bq, the peak is always located at 
negative Ba on the descending branch of the major loop. 
Such a peak position at negatiye Ba is a typical feature 
also in the parallel geometry.ErO However, it contrasts the 
case of a thin stripJn perpendicular field, where it was 
shown analyticallyE3 that for any Jc (i3)-dependence the 
peak is located exactly at Ba = 0. 

In the Bean model, there is a simple relation between 
the critical current and the width AAf of the major mag- 
netization loop. 



3(i , , ^ 
Jc = — AAf. 
2R 



(13) 



The same expression is often used to determine Jc from 
experimental AM data even when the width of the ob- 
served loop is not constant. As discussed in Refs. |^-^ 
the applicability range of such a procedure is limited. In 
the parallel geometry a simple proportionality only ap- 
plies for Ba larger than the full penetration field. For 
the thin disk case the field range where Jc oc AAf can be 
estimated from our calculations. Figure |^ shows Jc{B) 
inferred from the magnetization loop using Eq. (hq), to- 
gether with the actual Jc{B). One can see that at fields 
larger than the characteristic field 



Be — IJ'oJco/'^ 



(14) 



there is essentially no distinction between the two curves. 
We find that this holds independently of Bq and also for 
other Jc{B) models. Therefore, also for the present ge- 
ometry the i?-dependence of Jc can be inferred directly 
from AM{Ba), except in the low-field region. Here the 
correct Jc{B) can be obtained only by a global fit of the 
magnetization curve. 

The Bean-model virgin, nnagnetization for a thin disk 
can be expanded in Ba ast3E3 



fioM ~ XoBa 




(15) 



where xo = SR/^ird is the Meissner state susceptibility. 
Our numerical calculations show that the same expan- 
sion also holds for i?-dependent Jc, only with an effective 
value B^^ satisfying 



Bf/Bc « 1 - a ^/Bc/Bo. 



(16) 



We find that if Bq/Bc > 0.5 the parameter a = 0.50 
for the exponential model, and a = 0.43 for the Kim 
model. In the parallel geometry the low-field expansion 
has an additional B^ term which is not affected by the 
Jc{B) dependenceu Thus, the deviation from the Meiss- 
ner response at small Ba is there insensitive to Jc{B). 



This result contrasts the case of perpendicular geometry 
where due to demagnetization effects, a B-dependence of 
Jc affects the flux behavior even in the limit of low fields, 
see discussion in Ref. Ell. 



AM * 3d /2R 
actual J„ 




FIG. 3. The critical current inferred from the width, AM, 
of the major magnetization loop using Eq. ([L3|). For compar- 
ison the plot also shows the actual Jc{B) used in the calcu- 
lation: the exponential model with Bq = Be (lower curves) 
and -Bo ~ 3Bc (upper curves). The agreement is excellent for 
fields larger than B^. 



IV. COMPLEX AC SUSCEPTIBILITY 

A. Basic expressions 

The hysteretic dependence of the magnetization M as 
the applied field Ba is cycled leads to ac losses. The 
energy dissipation per cycle of Ba is 



W 



cycle 



. , dM(t) 
at 



(17) 



per unit volume. According to the critical state model, 
M{t) follows Ba{t) adiabatically, i.e., M{t) = M[Ba{t)]. 
Thus, the losses are given by the area of the magnetiza- 
tion hysteresis loop, § M dBa- 

It is conventional to express the ac response through 
the imaginary and real pacts of the so-called nonlin- 
ear magnetic susceptibility.E3 If the applied field is os- 
cillated harmonically with amplitude Barm i-6-, Ba{t) — 
Bam COS Lot, the magnetization is also oscillating with the 
same period. The complex susceptibility is then defined 
by the coefficients of the Fourier series of the in general 
anharmonic M{t), where the real and imaginary parts 
are given by 



Kn 



ttB, 



am Jo 



2Tr/LL 



M{t) cos{nu!t) dt 



ttB, 



am Jo 



2ixlu. 



M(t) sin(n,wi) dt , 



respectively. 

The dissipated energy, W , is determined by the re- 
sponse x^ at the fundamental frequency, namely 



x"^xi' 






2^0 



M{Ba)dBa. (18) 



Below we shall also analyze the real part of the suscepti- 
bility at the fundamental frequency, x' = x'j^, which can 
be expressed as 



X = 



2mo 



-Ba 



M{Ba) Ba dBa 
VBln-Bl 



(19) 



The x"{Bam) and x'{Bam) are calculated from these ex- 
pressions using M{Ba) obtained by the previously de- 
scribed procedure with Bam covering a wide range of 
amplitudes. For convenience, we normalize the suscepti- 
bilities to the Meissner state value xo — 8i?/37r(i.Ell 
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FIG. 4. Real (bottom) and imaginary (top) parts of the 
nonlinear susceptibility for a thin disk as functions of the am- 
plitude Bam of the applied ac field. Calculations are based on 
Eqs. (Ilq ) and (h9|) with Jc{B) given by the exponential model, 
Eq. (112), with Bo ~ SBc (curve 1) and Bq = Be (curve 2). 
For comparison the results for the Bean model are also shown 
(dashed line). 

As seen from Fig. ^ the response x" shows a maximum 
as a function of the field amplitude. Such a maximum 
is in fact a common feature in all geometries. For the 
Bean model for a long cylinder the peak is known to 
occur when Bam is equal to the full penetration field. 
In the perpendicular geometry the interpretation of the 
peak position is not so simple. Even in the Beanjnodel 
for a thin disk only numerical results are availabletil: the 
peak value equals Xmax ~ O-'iA and occurs at an ampli- 
tude of Bam — 1.94i?c, corresponding to the penetra- 
tion 1 — a„jR = 72%. We find that the B-dependence 



of Jc leads to a slight increase both in «„ and in the 
peak magnitude. For example, the numerical results for 
the Kim model with Be = Bq give Xmax — 0-29 and 
1 — am/R = 70%. The difference between various Jc{B) 
models becomes more distinct if one analyses the asymp- 
totic behavior at small and large field amplitudes, as 
shown below. 



B. Low-field behavior 

At small field -amplitudes the Bean model gives the 
exact expressionsOilj 



X'/Xo = -1 + lHBam/B,)y32, 

x"/xo - {Bam/B.f/n . 



(20) 
(21) 



Shown in Fig. a are our numerical results for x" for the 
exponential model. From the log-log plot it is clear that 
the quadratic dependence on Bam retains as in Eq. ( |2f| ) 
only with a modified coefficient. Moreover, we find that 
also V can be described by the Bean model expression 
Eq. ( po| ) with the same effective Be- The effective Be fits 
the expression 



Bf/B,^l-aBjBQ, 



(22) 



when Bq/Bc > 1 with a = 0.42 for the exponential 
model, and a — 0.36 for the Kim model. Interestingly, 
the same effective-description was found for the flux pen- 
etration depth af3 whereas it deviates from the descrip- 
tion of the virgin magnetization, Eq. ( |f q ) . 
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FIG. 5. Imaginary part of the susceptibility for the expo- 
nential model with Bq/Bc ~ 3 (curve f) and with Bq/ B^ ~ 1 
(curve 2). At low fields both curves as well as the dashed curve 
presenting the Bean model result follow the same quadratic 
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C. High- field behavior 

The high-field behavior of the dissipated energy W is 
shown in Fig. ^ for a variety of Jc{B) dependences. We 
choose to plot W rather than x" because the difference 
between the asymptotic behavior in the various models 



becomes more evident. One sees from the figure that 
for large Bam the Bean model yields W ex Bam- The 
exponential model shows saturation, whereas one finds 
after a closer inspection that the Kim model leads to 
a logarithmic increase. These behaviors can be under- 
stood by considering the fact that for large amplitudes 
the disk is fully penetrated and B{r) « Ba- Therefore, 
M{Ba) ex JciBa), and one obtains W ex J^"'" Jc{B) dB. 
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FIG. 6. High-field behavior of the dissipated energy, W ^ 
for the Bean, Kim and exponential models. 

The high-field behavior of the real part of the suscep- 
tibility, x', for different JciB) is shown in Fig. 0. For 
the Bean model we find asymptotically that x'/xo — 
— I.33(i3am/-Bc)~^^^ (dotted line), which is in agreement 
with Eq. 32 in Ref. O. For the _B-dependent Jc's we 
also find power-law behavior, although with different ex- 
ponents. For both the Kim and exponential model the 
asymptotic behavior is described by x' oc B/^/^. How- 
ever, also intermediate values for the exponent are pos- 
sible, e.g., for Jc = Jco/[l + {\B\/?>Bcfl'^\ the numerical 
results suggest that x' oc BaJ, ■ 

In order to understand this power-law behavior let us 
rewrite Eq. (n9[) as 



7i"S2,„ 



^oMrcv{Ba) Ba dBa 



vw. 



B^„ 



(23) 



where Mj-ov = Aff + M^ is the reversible magnetization. 
The integrand has different estimates in the regions I, II, 
and III indicated in Fig. ^. Therefore we divide the inter- 
val of integration correspondingly, x' — x'l + x'n + x'm- 
In region I, Mj-cv does not depend on Bam-, thus, x'l ^ 
B~^ at large Bam- In region II {Ba ^ Be) we use that 



M,ev(Sa) « / drr"" [J,{Ba + B,{r)) - J,{Ba - B,{r))] 



where Bi is the field created by the current. Expand- 
ing this expression one has Mj-ev oc J^{Ba) Jdrr^Bi{r). 



Then using the further simpHfication that Jdrr'^Bi{r) ex 
Jc{Ba), one obtains 



X7/OC 



^am J II 



J,{Ba) J',{Ba) Ba 



VW. 



52 



dB„ 



Taking Jc{B) oc (Bq/B)'' at large B, we arrive at the es- 
timates, x'li oc {Bc/Bamf (So/Bq„i)2s f^j. g^^g^n g^ and 
X'li oc (Be/Barn)^ (Bq/Bc)'^'' for large s. Finally, consider 
the region III, where Bam—Ba is of the order of /io>/c(Sa). 
Since the initial slope of the return branch does not de- 
pend on Bam, we have that M,-av{Ba) c>c Bam — Ba- It 

then follows that, x'lii oc {Bd Bamf {Bq/ Bamf ■ 

As the asymptotic behavior at large Bam is determined 

by the slowest decaying term, we arrive at the following 

result 

n-3(l+s)/2 -, 

^am , i ^ ± 

X « <; (24) 



It should be specially emphasized that the presented 
analysis for the high-field asymptotic behavior is not re- 
stricted to a thin disk. In fact, we expect the result 
( P4| ) to be valid in any geometry. This result is also in 
agreement with numerical calculations for long samples 
described by the Kim and the exponential modell 
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FIG. 8. Division of a magnetization loop into 3 regions 
treated differently when estimating x' a-t large Bam- 
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FIG. 7. High-field behavior of the real part of the sus- 
ceptibility, x') for the Bean model (curve 1) and for different 
Jc(B) dependences; Jc = Jco/(l + (|B|/3-Bc)^^^) (curve 2); 
Kim model (curve 3) and exponential model (curve 4), both 
with Bo — 3Bc- 

These power-laws fully agree with our numerical cal- 
culations shown in Fig. |^. The expression ^ gives the 
exact values for the exponent found for the Bean model 
(s = 0), the Kim (s = 1) and exponential (s = 00) mod- 
els and even for the Jc{B) with s = 1/2. Note however, 
that this asymptotic behavior is sometimes established 
only at rather low values of \x'\, see curve 4 in Fig. |^. 
Therefore one should be very careful in interpretation of 
corresponding experimental log-log plots. 



D. Plots of x" versus x' 



In contrast to graphs of x as a function of the field am- 
plitude or temperature, a plot of x" versus x' contains 
only dimensionless quantities, aiid||Ls,therefore very useful 
for analyzing experimental datar tM. In practice, such a 
parametric plot x"(x') can be obtained by scans either 
over the magnetic field amplitude or over the tempera- 
ture. Figure II presents the x"ix') plof of the data shown 
in Fig. 0. We observe that a i?-dependence of Jc gives 
a significant distortion of the graph. Compared to the 
Bean model one finds that: (i) the maximum is shifted 
to higher values of x" ', (ii) it occurs at smaller values of 
— x'; (iii) in the limit of large Bam (or high temperatures) 
x" falls to zero more abruptly. 
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FIG. 9. Parametric plot of the complex susceptibility for 
a thin disk using the three sets of curves presented in Fig. ]A. 
The dotted line shows the Bean model asymptotic behavior, 
Eq. (Eq). A B-dependence of Jc substantially distorts the 
plot, although not the slope at x' ^^ ~1- 



Meanwhile, at small Bam, as x' ~^ "1; the slope of 
x"ix') curve remains the same as in the Bean model, 
namely, 






32 

15^ 



1 



X_ 
Xo 



at x' 



-1 



(25) 



This result holds for any Jc{B). It also follows from the 
previous analysis showing that at low fields both x' and 
x" are modified by JAB) in the same way. The universal 
slope given by Eq. (|25| ) allows one to examine if exper- 
imental data are described by the critical state model 
without a priori knowledge of the actual Jc{B) depen- 
dence for the sample. 

The presented x"ix') plots for a disk in a perpendicu- 
lar field should be compared to similar plots for the long 
samples in a parallel field studied systematically in Ref. H. 
As expected, the Bean-model curve for a thin disk shown 
by the dashed line in our Fig. 9 appears quite differ- 
ent from the Bean-model curves for long samples shown 
in Fig. 7(a,b) of Ref. 0. Meanwhile, further analysis of 
these figures shows that the account of a B-dependent 
Jc always leads to very similar distortions of the x"(x') 
plots. Namely, in all geometries the x" peak increases 
in magnitude and shifts towards x' = 0. Note that such 
a behavior is found when the characteristic field Bq of 
the Jc (i?)-dependence is larger or of the order of Be- For 
Bf) ^ Be this behavior may change qualitatively. In par- 
ticular, in the parallel geometry, the peak nosition, Xmaxj 
becomes a nonmonotonous function of Bo.ti However, the 
case of Bq ^ Be is not very realistic for a thin disk since 
Be is proportional to the sample thickness while Bq is 
usually taken as geometry-independent. 
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FIG. 10. The behavior of x"(x') for various models. The 
Bean model predicts a peak located at Xmax = —0.38. A 
B-dependence in Jc shifts the peak to the right and changes 
the behavior at x' ^^ (our results), while flux creep shifts it 
to the left and changes the behavior at x' ^^ — 1 (Ref. |L9| ). 

It is interesting to compare our x"(x') plots to the ones 
obtained by calculations based on a non-linear current- 
voltage curve, j oc E^^", n < oo. Shown in Fig. p^ 
together with the CSM results is a x"(x')"Curve (dotted 



line) drawn in accordance to typical graphs presented in 
Refs. [lqjl9| . Compared to the Bean model curve, the 
maximum of x" increases in magnitude and shifts to- 
wards x' — —1- Moreover, the slope at x' — *■ — 1 becomes 
steeper. The last two features are in a strong contrast 
to the effect of having a S-dependent Jc in the CSM. 
Consequently, an analysis of the x"(x') plot allows one 
to discriminate between a strict CSM behavior and one 
where flux creep is an ingredient. 

Finally, we compare in Fig. |ll| our theoretical results 
to available expccimcntal data on the susceptibility of 
YBaCuO films O'Ej-lJ The shown data were obtained by 
reading selected points in the graphs found in the liter- 
ature. It is evident that the poor fit by the Bean model 
(dashed curve) is greatly improved by the curve (full line) 
calculated for a i3-dependent Jc- Whereas the agreement 
is better throughout the x"(x') plot; it is especially ev- 
ident at small |x'| (large field amplitudes), where the 
Jc (i3)-dependence plays a major role. There is still a 
discrepancy in the low-field region, where all experimen- 
tal points do not follow the universal CSM slope given by 
Eq. 25, The deviation can be caused by a flux creep lead- 
ing to a steeper slope.ESilj This suggestion can be checked 
experimentally by analyzing x"(x') plots obtained at dif- 
ferent temperatures. 







X'/Xo 
FIG. 11. Experimental susceptibilit*!-|data from Ref. |l^ 
(open circles) and Ref. |l^ (solid circles)Ej together with the 
CSM predictions for a thin disk: the Bean model (dashed 
curve) and the Kim model with Bo = 3Bc (full line). The 
Kim model gives a better agreement with experiment over 
the whole range. 



V. CONCLUSION 

Magnetization and ac susceptibility of a thin super- 
conducting disk placed in a perpendicular magnetic field 
were analyzed in the framework of the critical state model 
where Jc depends on the local flux density. We solved nu- 
merically the set of coupled integral equations for the flux 
and current distributions, and from that calculated mag- 
netization hysteresis loops as well as the susceptibility. 



X = x' + ix" ■ The results, which were obtained for sev- 
eral commonly used Jc decreasing with \B\, allowed us 
to determine the range of fields where the vertical width 
of the major magnetization loop, AAI{Ba), is directly 
related to Jc{Ba). 

We have shown that at small fields the virgin mag- 
netization and complex susceptibility have the same de- 
pendence on Ba as for the Bean model, although with 
difi'erent coefficients. For large ac amplitudes. Bam, the 
behavior of the ac susceptibility changes from x' oc Bam 
and x" oc B~^ for the Bean model, to x' '^ B~^ and 
x" oc Bam foi' Jc decreasing with \B\ as \B\~^ or faster. 
We could show numerically, and also presented an argu- 
ment, that when asymptotically Jc ~ |i?|~'', s < 1, one 

has x' Of Bam ^ ■ The results for the high-field be- 
havior of the susceptibility are expected to be valid for 
superconductors of any geometry. 

A most convenient test for critical-state models is pro- 
vided by an analysis of the x"{x') plot. We conclude 
that the asymptotic behavior at x' ^ — 1 is universal for 
the CSM with any Jc{B). On the other hand, flux creep 
can affect this behavior. The peak in x" at x' = —0.38 
predicted by the Bean model was found to be shifted to- 
ward x' = due to the _B— dependence in J^ and toward 
x' = — 1 because of flux creep. 
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